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Abstract 

The one dimensional dilute Kondo lattice model is investigated by means of bosonization for 
different dilution patterns of the array of impurity spins. The physical picture is very differ- 
ent if a commensurate or incommensurate doping of the impurity spins is considered. For the 
commensurate case, the obtained phase diagram is verified using a non-Abelian density-matrix 
renormalization-group algorithm. The paramagnetic phase widens at the expense of the ferromag- 
netic phase as the /-spins are diluted. For the incommensurate case, antiferromagnetism is found 
at low doping, which distinguishes the dilute Kondo lattice model from the standard Kondo lattice 
model. 

PACS numbers: PACS No. 75.10.Jm, 75.40.Mg, 05.50.+q 
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Heavy fermion systems have been of great theoretical interest since their discovery some 
twenty years ago.— The central problem posed by heavy fermion materials is to understand 
the interaction between an array of localized moments (generally / - electrons in lanthanide 
or actinide ions) and conduction electrons (generally s - or d - band). This situation is well 
described by an antiferromagnetically coupled Kondo type model. 

The solution of Kondo type models is well understood in two limiting cases; the single 
impurity limilj^ which can be reduced to a one dimensional problem and solved via Bethe 
Ansatz, and secondly the Kondo lattice model (KLM), which was solved via bosonization^ 
and numerous numerical approaches'"^ in one dimension for half - and partial conduction 
band - filling. For half - filling the results indicate the existence of a finite spin and charge 
gap. Accordingly in this case the Kondo lattice model is an insulator with well defined 
massive solitonic excitations of the spin sector. 

For partial conduction band filling, the conduction electrons form a Luttinger liquid, with 
spin and charge separation.^^ The localized spins, however, exhibit ferromagnetism, due to 
an effective double exchange coupling.^^ The double exchange is driving the system toward 
ferromagnetism, while the fluctuations generated by Kondo singlets compete against this 
tendency. As a consequence, the paramagnetic to ferromagnetic phase transition is of the 
quantum order - disorder type, typical to models with an effective random fieldjA However, 
for small Kondo coupling and close to half filling an RKKY liquid state and polaronic regime 
are always present.— 

Beyond these two solvable limits, no rigorous results exist for the intermediate cases, 
where the number of impurities are neither one, nor equal to the number of sites. This is the 
focus of our study. We concentrate on the one dimensional case, and start from the Kondo 
lattice limit introducing impurity spin holes, that is we will be dealing with a dilute Kondo 
lattice model (DKLM): 



where L is the number of sites and t > is the conduction electron hopping. We measure 
the Kondo coupling J in units of the hopping t. We denote by Nf {uf = Nf/L) the number 
(concentration) of impurities and A^^^ {uc = N^/L) the number (concentration) of conduction 
electrons. The constraint Nf < L is imposed by V, which is an operator that projects out a 
pre-determined set of /-spins. Sj are spin 1/2 operators for the localized spins, e.g. /, and 
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Si = 2^(^y ^i,cr^cry^i,cr' with CT the Pauli spin matrices and c-^., c^^^ the electron creation 
and annihilation site operators. 

We investigate the behavior of the DKLM both by an analytical approach, based on a 
standard bosonization scheme, and by numerical calculations. The latter were performed 
using the newly developed non-Abelian density-matrix renormalization-group (DMRG) 
algorithm,— which preserves the total spin and pseudospin symmetry. This choice of ba- 
sis first of all greatly facilitates the observation of magnetic phases and secondly it gives a 
dramatic performance improvement compared to the standard DMRG basis. 

The bosonization we use takes the standard approach^ by first decomposing the on-site 
operators into Dirac fields, with spinor components r = ± (otherwise known as the right, 
r = +, and left, r = — , movers): c^o- ~ J^r^^rxa = J2t ^^^''^'^T,a{x), where kp = nnc/2 
and we consider the lattice spacing to be unity. Next we bosonize the Dirac fields with 
\E't-_o- = exp(z$^o-)/v27rA, where 1/A is the ultraviolet cutoff. For the scalar Bose fields, 
^r,a{x) and its canonical conjugate momenta, Ur^^^x), ^r,a{x) = JUi^dx'Ilr^aix'), we use 
the standard Mandelstam representation,^^ which introduces a momentum cutoff function 
A(A;) = exp(— A|A;|/2) via the Fourier transforms. Thus, the electron field can be represented 
in terms of collective density operators which satisfy Bose commutation relations: 

c^,,^^ ^ exp{iTkFx) expi{9p{x) + r0p(x) + a[9^{x) + T(p^{x)]}/2 , (2) 

where the Bose fields (for both u = p,a) are defined by 4>u/0u = i{7i/N)J2kjto^^''^[^+{^^) =t 
i'_{k)]A{k)/k, where are the number fields and 61, the current fields. The charge (holon) 
and spin (spinon) number fluctuations are defined as Pr(^) = J2a PT,a{k), and ar{k) = 
J2a'^Pr,a{k)- Thls type of Bose representation provides a non-perturbative description of 
the conduction electrons in terms of® holons and spinous. We will neglect for the moment 
all the rapidly oscillating (umklapp) terms. These will give a contribution only at half filling, 
i.e. ric = rif, and will be analyzed later on. Thus, the bosonized form of DKLM is: 

+ ^ E {cos[0.(j)] + cos[2kFj + Mm {e~''^'-'^S+ + h.c.) 
j 

- A E MMj)] sin[2fc^j + ■ (3) 

2ttX ^ 

This equation has the same form as for a standard Kondo lattice^ii except i) we have to 
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keep in mind that the impurity spin, i.e. terms containing Sj, Sj' and 5"", contributes only 
if there is an / spin at site j, and ii) the even cutoff function A(fc), defined in Eq. Q, 
satisfying A(A;) ^ 1 for \k\ < 1/A and A(A;) ^ otherwise, is needed in the Bose fields to 
ensure that delocalized conduction electrons are described. Delocalization is essential to 
describe ferromagnetism^'^'^ since ferromagnetism in the Kondo lattice models is due to the 
double exchange, which only requires that Nc < A^/-— 

In this situation each electron has on average more than one localized spin to screen, 
and since hopping between localized spins is energetically most favorable for electrons which 
preserve their spin as they hop, this tends to align the underlying localized spins.— This also 
means that double exchange will vanish if the distance between impurity spins is larger than 
A. At lengths beyond A the electrons will behave as collective density fluctuations, as usual 
in one dimensional systems.— Hence, A measures the effective range of the double exchange, 
and in principle it is a function of J, ric and Uf. 

The most straightforward method to determine an ordering of the localized spins is by 
applying a unitary transformation H = e^He~^. We choose the transformation which 
changes to a basis of states in which the conduction electron spin degrees of freedom are 
coupled directly to the localized spins: S = ig^^ Z]j ^o-O) Sj. We perform the unitary 
transformation up to infinite order, so there is not any artificial truncation error generated 
(for details see Ref. O). In the new transformed basis the double exchange interaction 
leading to ferromagnetism is clearly exhibited and we obtain the effective Hamiltonian for 
the localized spins: 

J2 \ 

TT _ _ \ " qZ qZ 

+ 7^T.{^os[K{^)]+cos[2kF^]}S: 
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The K{i) is a term introduced by the unitary transformation, which counts all the Sfs to the 
right of the site i and subtracts from those to the left of i: K{i) = {J/2vp) J2iZi{S-_^i — S-_i). 
This term gives the crucial difference between KLM and DKLM, as will be explained later on. 
The most important term in Eq. is the first one, which clearly shows that a ferromagnetic 
coupling emerges for DKLM. This coupling is non-negligible for A^^^ < Nf and i ~ j < X and 
its strength will decrease with increasing distance between impurity spins. 



For Nc < Nf, the physical picture given by Eq. (@]) will be crucially different if the lattice 
of impurity spins contains commensurate or incommensurate array of holes. Hence, we 
analyze these two cases separately. If we have a commensurate doping of the impurity spins, 
then we can approximate the ferromagnetic term in the usual way by taking ^ l/uf for the 
shortest average distance between / spins: { J^?2jA/[27r^t>i?(H-njA^)]} J2i ^i^i+i/uf- Lattice 
sites which are not occupied by / spins are inert and do not contribute to the ferromagnetic 
phase. This was verified by DMRG: the calculated f — f spin correlation functions behave 
similarly as those of the normal KLM. The empty sites are inert. The /-structure factor 
has the usual peak at fc/vr = Nc/Nf for low J, hence in the commensurate case the DKLM 
behaves similarly to the standard KLM model. 

To understand the behavior of the second and third term from Eq. (@)), we notice that 
K{i) is vanishingly small for the commensurate the number of / impurity spins to 

the left and to the right of a given site i is the same. So the effective Hamiltonian will reduce 
to the random transverse field Ising model, as in the KLM. - The randomness is generated 
by (1 + cos[2kFi]) at large distances and it is driven by a cosine distribution, similarly to 
spin-glasses.— To determine the phase transition we need the dependence A = A( J, rzc, nj), 
which however is very difficult to determine and as such we use the low density value A ~ 



y'2/J close to criticality, similarly to previous works. ^'^^ In this way we obtain the critical 
phase transition line to be: J = tt sin(7™c/2)/[l — tt sin{7mc/2)/n'j], which represents a 
quantum order - disorder transition with variable exponents.^ However, this ferromagnetic 
phase disappears for larger distances between impurities because, as mentioned earlier, the 
double exchange interaction vanishes if the average distance between impurity spins, l/'n,/, 
is larger than A. This is very important because it ensures that the single impurity limit, 
rij — > 0, is free of ferromagnetism, as it should be. 

The incommensurate case is more difficult than the commensurate case. The reason is 
that in the low concentration limit the properties of DKLM will be very much dependent 
on the random distribution of / spins. We may observe phase separation or clusterization 
processes in this case. In this limit, where the average distance between impurities is very 
large, then the single impurity^ approximation seems natural. However, if we look at small 
doping of / electrons only, then the main difference compared to the commensurate limit 
studied previously is that the K{i) term, in Eq. Q, is not negligible anymore. The impu- 
rity / spins are no longer equally distributed to the left and right of a given site i, hence 
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K{i) ^ (—!)'( J/2t'ir), which gives rise to a staggered field. The properties of Eq. (jH) are then 
given by the staggered field Ising model, which gives an antiferromagnetic ordering. This 
antiferromagnetic ordering of the impurity spins represents a new element in DKLM com- 
pared to the Kondo lattice. This corresponds to the soliton lattice obtained by Schlottmann 
in a dynamical mean-field treatment of the three dimensional dilute Kondo lattice.^^ 

Similar behavior also occurs above half- filling, i.e. Nc > Nf, where double exchange (as 
shown previously) does not appear. But bosonization still works: the effective Hamiltonian 
reduces to the second and third terms of Eq. from which the most dominant term, for 
low doping of impurity spins, as in the case described previously, is a staggered S'f field. 
As the first term in Eq. (@)) is missing in this case, the only fiuctuation which can destroy 
a locked staggered order is S^. For large J (4 < J) the staggered order wins. While for 
smaller values of J the systems will be disordered. 

As we approach half filling from both sides, the bosonization approach breaks down 
as the strongly oscillating (umklapp) fields start dominating. The DKLM will undergo 
a metal-insulator transition as in a standard quantum sine-Gordon model^^ by dynamical 
mass generation. A Spin gap will also appear. This can be understood easily, because the 
half filled DKLM is equivalent to the quarter filled periodic Anderson model, which has a 
antiferromagnetic order.>i2. The only difference from the Kondo lattice is that the massive 
solitons obtained for DKLM are of Su, Schrieffer and Heeger type.— 

To confirm the previously obtained magnetic phases, we performed non-Abelian DMRG 
analysis of the DKLM model in the commensurate case, for both the Nc < Nf and Nc > Nf, 
where ferro- and antiferromagnetism, respectively exist. For the Nc < Nf case, we have 
investigated several commensurate dilution patterns of the form: "0 / / /..."; 
"0 / / / / . . ."; and "0 // // //..." in a 64-site long chain. While in 
the Nc > Nf limit we have studied the n/ = 0.8 localized spin filling (to be in accordance 
with the bosonization requirement of low dilution) on a 80-site long DKLM chain, i.e., we 
investigated the pattern: "0 / / / / / / / / / . . .". 

These patterns were selected in such a way that the chain middle-point refiection sym- 
metry was preserved to accelerate the calculations. There was only one exception: the "0 
/ / /" chain has an impurity in the middle: "0 / / / / 0" , but its effect is rather 
small compared to our final errors and so it was neglected. 

In the DMRG calculations careful error and convergence analysis were used, and we 
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extrapolated the energy linearly to zero truncation error (we saw no quadratic terms large 
enough to affect the fit). For each given dilution pattern, filling Nc/Nf, and interaction 
constant J, we used several DMRG sweeps of between 200 and 400 S'0(4)-symmetric states. 

For Nc < Nf ferromagnetism appears at large J values, see Fig. ^ A point on the phase 
diagram shown in Fig. Q is judged to be ferromagnetic if the extrapolated energy of the 
spin S'max run is lower than the spin zero energy. This energy difference to the spin singlet 
excited state can be calculated directly using the 5*0(4) basis set of the non-Abelian DMRG. 
The phase transition line can be determined with high accuracy if one plots the energy gap 
between the spin S'max and spin zero states as a function of J - the gap rapidly decreases as 
the transition is approached. It should be mentioned that, as well as in the standard Kondo 
case,— there also exists additional ferromagnetic phases inside the paramagnetic region for 
dilute Kondo chains. 

In the opposite limit, i.e., > Nf, we have confirmed the existence of antiferromag- 
netism by calculating the spin structure factor S{k). As can be seen in Fig. |2l for small J 
and away from half filling the peak in S{k) is at 2/ci? (similarly to KLM^), i.e., the DKLM 
is a disordered paramagnet. The dominant 2k p backscattering processes are manifest of a 
system of free localized spins embedded in effective fields determined by conduction electron 
scattering. The half filled case, i.e., nc = 1 as in the KLM mode)^ is always an antifer- 
romagentic insulator. The tendency of Kondo lattice models to form an antiferromagnetic 
insulator at half filling for any nonzero J is a consequence of the Kondo effect.— As we 
increase J, the distinctive antiferromagnetic peak at k = n appears away from half filling 
(see Fig. 121). As we further increase J the 2kF peak broadens and decreases in strength while 
the antiferromagentism becomes more robust, at J = 5 it appers already at ^ 0.925. 

In conclusion, we have studied the dilute Kondo lattice model in one dimension both 
numerically, using DMRG, and analytically, with a standard bosonization approach. We 
have derived an effective Hamiltonian for the / spins, which accounts for the appearance 
of a ferromagnetic phase seen with a commensurate dilution pattern of the impurity spin 
array. The paramagnetic-ferromagnetic phase transition shifts to higher coupling J values 
as the /-spins of the chain are diluted, in agreement with the numerical DMRG calculation. 
We have also shown that for incommensurate dilution or above half filling, i.e., A^^^ > Nf 
for low doping of impurity spins, antiferromagnetism is found. This distinguishes the dilute 
Kondo lattice model from the standard Kondo lattice model. 
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FIG. 1: The phase diagram of the dilute Kondo model in different commensurate filling cases 
for Nc < Nf. Legend shows patterns of dilution. Open circles correspond to the standard KLM 
model. The system of a given dilution pattern is ferromagnetic above the corresponding solid line 
and paramagnetic below. 
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FIG. 2: Typical J dependences of the conduction electron spin structure factor, S{k) for n/ = 0.8 
and ric = 1 continuous, Uc = 0.975 dotted, ric = 0.95 dashed, ric = 0.925 long dashed and ric = 0.90 
dot-dashed curves, respectively. 
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